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ABSTRACT 

We present the first results of an ongoing spectroscopic follow-up of close luminous red galaxy 
(LRGs) and Mgll AA 2796,2803 absorber pairs for an initial sample of 15 photometrically selected 
LRGs at physical projected separations p < 350 kpc from a QSO sightline. Our moderate- 
resolution spectra confirm a physical association between the cool gas (T ~ 10^ K) revealed by the 
presence of Mgll absorption features and the LRG halo in five cases. In addition, we report an 
empirical estimate of the maximum covering fraction (/^max) of cool gas in massive, > 10^^ Mq 
dark matter halos hosting LRGs at z ~ 0.5. This study is performed using a sample of foreground 
LRGs that are located at p < 400 kpc from a QSO sightline. The LRGs are selected to have a 
robust photometric redshift az/{l + z^h.) 0.03. We determine /^max based on the incidence of Mgll 
absorption systems that occur within Zph ± 3 cr^ in the spectra of the background QSOs. Despite the 
large uncertainties in 2;ph, this experiment provides a conservative upper limit to the covering fraction 
of cool gas in the halos of LRGs. We find that A^max ^ 0.07 at Wr{2796) > 1.0 A and /^max ^ 0-18 at 
^^(2796) > 0.5 A, averaged over 400 h ^ kpc radius. Our study shows that while cool gas is present 
in > 10^ Mq halos, the mean covering fraction of strong absorbers is no more than 7%. 
Subject headings: Quasarsiabsorption lines — Cosmologyidark matter — galaxies:evolution 



1. INTRODUCTION 

A detailed description of how galaxies acquire their gas 
is essential to establish a comprehensive theory of galaxy 
evolution. Analytical calculations and hydrodynamical 
simulations have shown that the gas in M > 10^^ Mq 
halos^, is heated to high temperatur e by yirial shocks 



e.g., Birnboim & Dekel 2003; Keres et al. 2005 



Dekel 



|06 Birnboim||2006j |Keres et al.|,2009^ . Within the hot 
gas halo, thermal instabilities can induce the forma- 
tion of pressure-supported clouds of cool gas (e.g., Mo| 
[fc Miralda-Escud e 1996; Mailer fc Bullock 2004). Such 
clouds have been resolved In high-resolution hydrody- 
namic simulations of Milky- Way type dark matter ha- 
los (e.g., Kauffmann et al. 2008; Keres & Hernquist 
2009). These clouds bring some baryonic material into 
the galaxy, but the overall cooling rate is reduced be- 
cause the hot phase has lower density. Whether or not 
these clouds can form and survive in more massive halos 
is an open question. 

Our empirical knowledge of the cool gas content of dark 
matter halos has been shaped pri marily by HI 21-cm 



observations of lo cal galaxies (e.g., Thilker et al. 



2004 



Doyle et al.||2005[ |Verdes-Montenegro et al |[200Tfl2007p . 
At 2; > 0.1, we rely on intervenine; systems detected alone; 



At 2; > 0.1, we rely on intervening systems detected along 
the sight lines to background QSOs to describe the gas 
content of dark matter halos. (e.g., Lanzetta & Bowen 
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19901 ISteidel et al.|[T994l |Chen et al.|[2QQT| 
ker 2008). A number of studies have used 



^ ^ ^ . Chen fc Tin- 

ker 2008). A number of studies have used the Mgll AA 
z796,2803 absorption doublets to trace cool, T ~ 10^ K, 
gas and its physical association wi th foreground galaxies 



located near QSO sightlines (e.g., Tripp fc Bowen||2005^ 
Kacprzak et al.||2008[ |Chen fc Tmker 2008]). Insights 



into the cool gas content of dark matter halos can also 
be obtained by measuring the large-scale clus t ering am 
plitude of Mg II systems (Bouche et al. 2006 Tinker & 



Chen 



2008 



^[2QTq|. 

nier et 



Gauthier et al. (|2009[), we calculated the cluster- 



ing amplitude of strong Mg II absorbers (with rest-frame 
absorption equivalent width Wr(2796) > 1 A) using lu- 
minous red galaxies (LR Gs) at z ^ 0.5. LR Gs are 
old and passive galaxies ( jEisenstein et al.l |200ip resid- 
ing in M > 10^^ Mq halos. (e.g., Zheng et al 



2008 



Blake et al.| [20081 [Padmanabhan et al. 2008| [GautEier 
et al.||20Q9 ). They are identified using photometric red- 



shift techniques that offer a ty pical redshift a ccuracy of 
cr^ ^ 0.045 at Zph 0.5 (|Colli'ster et al. | |2007|. 

An interesting result trom [Gauthier et al. ( [2009[ ) is 
that the LRG-Mg II cross-correlation signal is compara- 
ble to the LRG auto-correlation on small scales (< 300 
CO- moving kpc ) that are well within the virial radii 
of the halos. The results suggest the presence of cool gas 
inside the dark matter halos of LRGs. 

Here we report the first results of an ongoing spec- 
troscop ic follow-up of the close LRG-Mg II pairs found 
in Gauthier et al. (2009). In five of the 15 cases, the 
precise spectroscopic redshifts establish a physical con- 
nection between these LRG-Mg II pairs. The spectro- 
scopic study of LRGs is supplemented with a survey of 
Mg II absorbers in the vicinity of photometrically identi- 
fied LRGs in the SDSS data archive. This survey allows 
us to utilize the vast survey data available in the SDSS 
archive to derive additional constraints on the covering 
fraction of cool gas in 10^^ Mq halos. We adopt a A 



2 



Gauthier et al. 



cosmology with Qm = 0.3 and ft a = 0.7, and a dimen- 
sionless Hubble parameter h = i^o/(100 km Mpc~^) 
throughout the paper. All distances are in physical units 
unless otherwise stated. 



2. EXPERIMENTS 

To examine the cool gas content in massive halos, 
we have designed two experiments. The first one is a 
spectroscopic follow-up study of photometrically identi- 
fied LRGs located near Mgll absorbers (Sample A be- 
low). The spectroscopic observations allow us to identify 
physically associated LRG-Mg II absorber pairs for inter- 
preting the small-scale clustering signal seen in | Gauthier] 
[eFal.l (|2QQ9|) . In addition, the optical spectra offer addi- 
tional knowledge for the stellar and ISM properties of the 
LRGs. The second experiment is a survey of Mgll ab- 
sorbers in the vicinity of photometrically identified LRGs 
(Sample B below). This study allows us to make use of 
the vast amount of galaxy and QSO survey data already 
available in the SDSS archive and derive statistically sig- 
nificant constraints for the covering fraction of cool gas 
around LRGs. 

2.1. Sample A - LRGs in the vicinity of known Mgll 
absorbers 

The purpose of the spectroscopic study of LRGs is 
to establish a physical association between the cool gas 
traced by the Mgll absorbers and the dark matter ha- 
los of the LRGs. The LRGs are selected based on the 
photometric redshifts and typical photometric redshift 
uncertaintie s for these galaxies i s = 0.03(1 + %hot) 
for i' = 19 (ICollister et al.||2007|). The identification of 



physical pairs requires precise spectroscopic redshifts for 
the galaxies. 

To achieve this goal, we established a catalog of close 
Mgll-LRG pa irs for follow-up spe ctroscopy by cross- 
correlating th e|Prochter et aL| ( 2QQ6[) Mgll absorber cat- 
alog with the 'Collister et al.^ ('2007") MegaZ-LRG sam- 
ple. The Mgll catalog is an extension of the SDSS DR3 
sample to include DR5 QSO spectra. This sample of 
Mgll absorbers has a 95% completeness for absorbers 
of Vl^^(2796) > 1 A . The catalog contains 11,254 ab- 
sorbers detected at ^^Mgii = 0.37 — 2.3 along 9,774 QSO 
sightlines. We excluded absorbers within 10,000 kms~^ 
from the QSO redshift to avoid associated absorbers of 
the QSO. We excluded absorbers with 2;Mgii < 0.4 or 
^Mgii > 0.7 to match the redshift interval of the LRGs. 
A total of 2461 Mgll absorbers satisfied our selection 
criteria. 

The MegaZ-LRG catalog is a photometric redshift cat- 
alog of approximately 10^ LRGs found in the SDSS DR4 
imaging footprint. The catalog covers more than 5,000 
deg^ in the redshift range 0.4 < z < 0.7. When cross- 
correlating the two catalogs, we rejected LRGs located 
within the photometric redshift uncertainty 3 x from 
the QSO redshift. We selected only those LRGs that are 
in the foreground of the QSO. 

We found a total of 646 Mgll-LRG pairs within pro- 
jected separation p = 400 (physical) kpc . Note that 
the sample of 64 close Mg II-LRG pairs with p < 300 
kpc found in Gauthier et al. (2009) is included in 



only kept isolated LRGs.^. Our final catalog for follow- 
up galaxy spectroscopy consisted of 331 'isolated' LRGs 
that are near a known Mgll absorber. We call this sam- 
ple Sample A. 

The primary utility of Sample A was to understand 
the sm all-scale clustering signal seen in Gauthier et al.| 
( 2009 ) . The sample of spectroscopically confirmed phys- 
ical LRG-Mg II pairs also allows us (1) to examine pos- 
sible correlations between the stellar properties of the 
LRGs and the presence of cool gas at large radii, and (2) 
to study the kinematics of cool gas in massive halos with 
follow-up echelle spectroscopy of the absorbers. 



2.2. Sample B - random LRGs near QSO sightlines 

To derive constraints on the covering fraction of cool 
gas, we first establ ished a QSO-LRG pa ir sample by 
cross-correlating the Schneider et al. (|2007[) QSO c atalog 
and the MegaZ L RG ca talog (jCollister et al.'20071). The 
[Schneider et all (|2007|) catalog contains more than 77K 
QSOs in the SDSS DR5 archive. MegaZ is a photometric 
redshift catalogue of approximately 10^ LRGs in the red- 
shift range 0.4 < Zph < 0.7. The redshift uncertainties 
are estimated to be o-z/{l + %h) = 0.03. 

We considered QSOs that are in the distant back- 
ground from the LRGs with (^qso — ^ph) > 3 in order 
to exclude correlated QSO-LRG pairs. In addition, we 
included only QSOs with ^Qgo < 18-5 and at z < 1.45, 
in order to obtain sensitive limits for Wr{2796) and to 
remove contaminations due to C IV absorption features. 
Furthermore, we excluded LRGs found in groups of two 
or more neighbors. Finally, we considered only pairs with 
p < 400 kpc at Zph, a maximum separation that is 
about the expected size of the host dark matter halos. 
This procedure yielded a total of 620 LRG-QSO pairs. 
We call this pair sample Sample B. This catalog contains 
no a priori information about the presence of Mg II ab- 
sorber in the QSO spectra. It was designed to compute 
the covering fraction of cool gas in the dark matter halos 
of LRGs (see the analysis in § 4.2). 



this larger sample of 646 pairs. From this sample, we 



3. OBSERVATIONS AND DATA REDUCTION 

In May 2009, we obtained medium-resolution spec- 
tra of three LRGs in Sample A (all within p = 160 
kpc of a known Mgll absorber) using the Double 
Imaging Spectrograph (DIS) on the 3.5 m telescope at 
the Apache Point Observatory. The blue and red cam- 
eras have a pixel scale of 0.4" and 0.42" per pixel, respec- 
tively, on the 3.5 m telescope. We used the B400/R300 
grating configuration with a 1.5" slit. The B400 grat- 
ing in the blue channel has a dispersion of 1.83 A per 
pixel and the R300 grating in the red channel has a dis- 
persion of 2.31 A per pixel. The blue and red channels 
with the medium resolution gratings together offer con- 
tiguous spectral coverage from A = 3800 A to A = 9800 
A and a spectral resolution of FWHM ^ 500 kms~^at 
A= 4400 A and FWHM ^ 400 kms"^at A = 7500 A 

^ An isolated LRG is one that does not have a a neighboring 
LRG within the volume defined by their virial radius and photo- 
metric redshift uncertainty. This approach allow us to reduce the 
ambiguity of attributing the Mgll absorber to one or more LRGs. 
We address possible contamination due to blue satellite galaxies in 
§5. The treatment of groups is beyond the scope of this letter and 
will be discussed in a forthcoming paper (Gauthier et al. 2010 in 
preparation) 
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.The observations were carried out in a series of two 
exposures of between 1200 s and 1800 s each, and no 
dither was apphed between individual exposures. Flat- 
field frames were taken at the end of the night for each 
night. Calibration frames for wavelength solutions were 
taken immediately after each science exposure using the 
truss lamps on the secondary cage. The DIS observa- 
tions are characterized by a typical seeing of 1.3" (as 
measured with the slitviewer guide star). Note that we 
also retrieved the spectrum of one LRG in Sample A 
(SDSSJ113731.00+060748.6) from the SDSS archive and 
included it in our spectroscopic studies. 

We pursued our spectroscopic observations of LRGs 
in September 2009 using the Boiler & Chivens (B&C) 
spectrograph on the du Pont 2.5-m telescope at Las 
Campanas Observatory, Chile. We obtained long-slit 
medium-resolution spectra for seven additional LRGs in 
Sample A. In addition, we selected four LRGs from sam- 
ple B for the spectroscopic B&C observations. The B&C 
camera has a plate scale of 0.7" per pixel. We used the 
300 In/mm grating with central wavelength of 6500 A 
and grating angle of 6 degrees. Since most spectral fea- 
tures of early- type galaxies at z ~ 0.5 are located at 
A > 5000 A , we employed the blocking filter GG4495 
to block the light blueward of ^ 5000 A , in order to 
avoid second-order contamination from the blue part of 
the spectrum. The slit width for all observations was 
1.1". The adopted grating configuration had a disper- 
sion of 3 A per pixel and a wavelength coverage of more 
than 6000 A . For the calibration frames, we followed the 
methodology employed for the DIS observations. A series 
of dome and sky fiats were taken each afternoon prior to 
the observations and calibration frames for wavelength 
solution were taken immediately after each science expo- 
sure. The typical seeing measured with the guide star 
was 1-1.4" during most observations. 

All of the DIS and B&C spectroscopic data were re- 
duced using standard long-slit spectral reduction pro- 
cedures. The spectra were calibrated to vacuum wave- 
lengths, corrected for the heliocentric motion, and fiux- 
calibrated using a spectrophotometric standard. Red- 
shifts of the LRGs were determined based on a cross- 
correlation analysis using known SDSS templates. The 
typical redshift uncertainty is ~ 0.0003. 

In summary, the spectroscopic sample consists of 11 
galaxies in Sample A and four in Sample B. A journal 
summarizing the observations of the 15 galaxies can be 
found in Table 1. 

4. RESULTS 
4.1. Optical spectra of LRGs 

The optical spectra described in § 3 confirmed that 
five of the LRGs occur within 350 kms~^ of a Mgll ab- 
sorber. Four of these LRGs belong to sample A and one 
to sample B. Figure 1 shows the optical spectra and SDSS 
images of these five LRGs. The close proximity of these 
LRGs and Mgll absorber pairs (p < 300 kpc and 
lA'^l < 350 kms~^) strongly argues for a physical asso- 
ciation between the galaxy and the absorbing gas. We 
note that the volume spanned by \Av\ = ±350 kms~^ 
and p = 350 kpc is ~ 5 co- moving {h~^ Mpc )^ 
at z = 0.5. Within this small volume, we estimate 
based on the best-fit luminosity functions determined 



for red galaxies at z ^ 
d)^ = 5.16 X 10-^ {h-^ 



- 0.5 (M* 
Mpc 



5 log h -19.8: 



see e.g., 



Brown et al. 



2007||2008P that the probability of the LRG being a ran- 
dom galaxy is negligible. 

In Figure [2] we present absorption profiles of the Mgll, 
Mgl and Fell transitions from the SDSS QSO spectra 
for these five physical pairs. In each panel, zero veloc- 
ity corresponds to the spectroscopic redshift of the LRG. 
All absorbers occur within \Av\ < 350 kms~^ of the 
spectroscopically identified LRGs. Although the spec- 
tral resolution of SDSS spectra does not allow us to re- 
solve individual components, the panels show that three 
of the Mg II absorbers have associated Fe II and Mg I ab- 
sorption features. 

All 15 LRG spectra exhibit spectral features domi- 
nated by absorption transitions, indicating an old un- 
derlying stellar population and little star formation in 
the recent pas t. A stellar popula tion synthesis analysis 
done with the Bruzual & Chariot] ([2003 ) spectral library 
established that the most likely stellar population ages 
for these LRGs span 1-11.75 Gyr (Gauthier et al. 2010 
in preparation). 

4.2. Incidence of cool gas in LRG halos 

We determined the covering fraction of cool gas in mas- 
sive halos, based on the presence/absence of Mgll ab- 
sorbers in the spectrum of the background QSO near the 
redshifts of the LRGs in Sample B. For each LRG-QSO 
pair, we visually inspected the QSO spectrum in search 
of Mgll absorbers in the redshift uncertainty range of 
the LRG, Zph ± 3az. For ^ 0.045 at Zph = 0.5, this 
interval corresponds to approximately ±400A centered 
at ~ 42OOA . We defined this interval as the search inter- 
val, which is large enough to identify all possible Mgll 
absorbers that may be associated with the LRG. De- 
spite the large redshift uncertainties, however, our search 
yielded a conservative maximum gas covering fraction 
i^max of Mgll absorbing gas in LRG halos. If no ab- 
sorber was visually detected, we measured a 2-a upper 
limit per spectral resolution element, Wiim, for the equiv- 
alent width of the Mgll A 2796 transition, using the noise 
spectrum. For each sightline with a detection, we mea- 
sured the redshift and strength of the absorber based 
on a Gaussian profile analysis. For these sightlines, we 
also determined Wum following the procedure described 
earlier. Because the S/N varies substantially between 
QSO spectra, VFum represents the minimum absorption 
strength that can be recovered in each QSO spectrum. 
It serves to define a homogeneous sample of QSO spectra 
that can be adopted to determine the incidence of Mg II 
absorbers. 

To determine the incidence of Mgll absorbers with 
Wr{2796) > Wo A around LRGs, we first identified from 
the parent LRG-QSO pair sample the QSO spectra that 
satisfy Wum < Wo A to ensure that the QSO spectra have 
sufficient S/N for revealing a Mgll absorber of minimum 
strength Wq. Next, we determined /^max by evaluating 
the fraction of the QSO spectra that display a Mgll ab- 
sorber of Wr {2796) > Wo in the vicinity of the LRGs. 
We calculate /^max for different values of Wq in different 
p intervals. 

This search allowed us to estimate the maximum pos- 
sible covering fraction of cool gas based on the incidence 
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TABLE 1 

Summary of the Long-slit Spectroscopic Observations of Luminous Red Galaxies 



ID 


RA(J2000) 


Dec(J2000) 


-^phot 




Instrument 


Exptime (sec) 


UT date 


Sample 


SDSSJ003816.29-092550.5 


00:38:16.29 


-09:25:50.5 


0.51 


18.70 


B&C 


2x1800 


2009- 


-09- 


21 


B 


SDSSJ011942. 14-090225.4 


01:19:42.14 


-09:02:25.4 


0.49 


19.88 


B&C 


3x2400 


2009- 


-09- 


18 


A 


SDSSJ015452.46-095533.6 


01:54:52.46 


-09:55:33.6 


0.55 


19.98 


B&C 


2x2400 


2009- 


-09- 


21 


B 


SDSSJ021819.25-083331.8 


02:18:19.25 


-08:33:31.8 


0.56 


19.46 


B&C 


2x2400 


2009- 


-09- 


19 


A 


SDSSJ023705.35-075513.7 


02:37:05.35 


-07:55:13.7 


0.55 


19.44 


B&C 


1800 + 900 


2009- 


-09- 


21 


B 


SDSSJ034802.50-070339.3 


03:48:02.50 


-07:03:39.3 


0.49 


19.08 


B&C 


2x2400 


2009- 


-09- 


18 


A 


SDSSJ113731. 00+060748.6 


11:37:31.00 


+06:07:48.6 


0.61 


18.00 


SDSS 










A 


SDSSJ142610.27+594704.6 


14:26:10.27 


+59:47:04.6 


0.48 


18.69 


DIS 


1800 + 1500 


2009- 


-05- 


26 


A 


SDSSJ160725.87+471221.7 


16:07:25.87 


+47:12:21.7 


0.56 


19.70 


DIS 


2x1800 


2009- 


-05- 


26 


A 


SDSSJ161713.60+243254.8 


16:17:13.60 


+24:32:54.8 


0.57 


19.06 


DIS 


2x1800 


2009- 


-05- 


26 


A 


SDSSJ204820.88-001640.1 


20:48:20.88 


-00:16:40.1 


0.45 


19.47 


B&C 


2x2400 + 1800 


2009- 


-09- 


19 


A 


SDSSJ205145.36-052121.2 


20:51:45.36 


-05:21:21.2 


0.50 


19.44 


B&C 


2x2400 


2009- 


-09- 


18 


A 


SDSSJ211625.92-062415.4 


21:16:25.92 


-06:24:15.4 


0.55 


19.31 


B&C 


2x2400 + 2000 


2009- 


-09- 


21 


B 


SDSSJ212713. 92-000747.2 


21:27:13.92 


-00:07:47.2 


0.47 


19.19 


B&C 


2x1800 


2009- 


-09- 


20 


A 


SDSSJ232450.16-095048.4 


23:24:50.16 


-09:50:48.4 


0.53 


19.79 


B&C 


2x2400 


2009- 


-09- 


19 


A 




1 

0.5 




= 1 1 = : 4, 

- « ; X 

— o: : 


1 o 


1 1 1 I.I 1 loo 


1 ' 


SDSSJ2I 1625.92-06241 5.4 

20" 

,/ ■ 

Wr(2796)=0.63 A Z=0.5240 
|Av| = 79km/s lOOkpc/h 




^ l' 1 !: 1 


1 


1 1 1 1 1 1 III 1 1 - 



6000 7000 

Wavelength (A) 



8000 



Fig. 1. — Spectra of the LRGs in our spectroscopic sample with associated Mgll absorbers, along with the thumbnail images of the 
galaxies. The top four LRGs were selected from Sample A (§ 2.1) and the bottom one was selected from Sample B (§ 2.2). For each LRG, 
we show the reduced spectrum in thick solid histograms and the corresponding 1-a error array in thin solid histograms. The dotted features 
are contaminating sky lines or artifacts. The thumbnail images are reproduced from the SDSS data archive to show the relative alignment 
of the LRG-QSO pairs. The LRG is at the center of each image and is connected to the background QSO by the arrow. The redshift of 
the LRG and the projected distance to the QSO sightline are listed in the bottom right corner. The velocity separation between the LRG 
and Mgll absorber is listed in the bottom left along with Wr (2796). 
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Sample A 



Sample B 



I ZLRG=0.5409 ZLRG=0.4976 zlrg=0.5703 zlrg=0.5446 | 

SDSS J1 42609.74+594629.5 SDSS J 160726.77+471 251 .3 SDSS J1 61 71 4.1 2+243255.6 SDSS J232448.00-095142.3 




-500 500 -500 500 -500 500 -500 500 

Relative Velocity (km/s) 



I zlrg=0.5240 

SDSS J21 1626.40-062437.4 



Mgn 27dB 




Mgn 2803 



1 



Mgl 285^. 




If b 



Fen 2100- 



-500 500 



Fig. 2. — Absorption profiles of the five physical LRG-Mgll absorber pairs in our spectroscopic sample of LRGs. The first four objects 
belong to Sample A and the last one to Sample B. We included the associated Mgl A 2852 and Fell A 2600 transitions. Zero velocity 
corresponds to the spectroscopic redshift of the associated LRG. Contaminating features are shown in gray. 
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Fig. 3. — The incidence of cool gas as probed by the pres- 
ence/absence of Mgll absorption features, ^Cmax, at different pro- 
jected radius of massive LRG halos in Sample B. The vertical 
dashed line represents the typical virial radius of the LRG pop- 
ulation (see Gauthier et al. 2009). Despite the large redshift un- 
certainties of the LKGs, the observed incidence of Mgll absorbers 
places a strong constraint on the maximum covering fraction of cool 
gas around LRGs. We estimate the contribution due to random 
background absorbers that occur within the redshift interval of the 
LRGs by chance coincidence. The result is shown in the dotted line. 
Including correlated absorbers due to large-scale matter cluster- 
ing does not increase the background contribution to beyond 3%. 
The solid curve corresponds to the expected maximum contribution 
from satellite galaxies to the observed incidence of Wr (2796) > 1.0 
A absorbers (see Discussion section). Our study shows that nei- 
ther satellite galaxies nor background random/correlated galaxies 
are sufficient to account for the observed incidence of Mgll ab- 
sorbers in the vicinity of LRGs. Note that the errobars on ^tmax 
are derived from the Poisson distribution. 

of Mgll in LRG halos (see § 5 for a discussion on the 
contamination rate due to correlated and random galax- 
ies). We first considered the incidence of Mgll absorbers 
at different absorption threshold Wo- All QSO spectra 
in the 620 close LRG-QSO pairs have sufficient S/N 
{W\im < 1 A) for uncovering an absorber of WV(2796) > 
1 A. Our search showed that 45 of the 620 LRGs have 
a Wr{2796) > 1 A absorber at p < 400 h'^ kpc, indi- 
cating a mean covering fraction averaged over the entire 
LRG halos of /^max(p < 400) = 0.07 for W^(2796) > 1 A 
absorbers. Next, we calculated /^max in different p inter- 
vals. Figure [s] shows that for Wr ^ lA Kma^{p) ^ 10% 
at p < 300 h kpc but declines to /^max ^ 5% at larger 
p (triangles). 

Similarly, 575 QSO spectra in the 620 close LRG-QSO 
pairs have sufficient S/N {Wum < 0.5 A) for uncovering 
an absorber of W^(2796) > 0.5 A, and 105 of the 575 
LRG members have a Wr(2796) > 0.5 A absorber at 
p < 400 kpc. The search therefore yielded a mean 
covering fraction averaged over the entire LRG halos of 
/^max(/> < 400) = 0.18 for absorbers of W^(2786) > 0.5 
A. The estimated /^max versus p for Wr(2796) > 0.5 A 
absorbers are also presented in Figure [s] (squares). 

The ^max results at small separations [p ^ 100 
kpc ) are consistent with the covering fraction mea- 



surements from the spectroscopic Sample B. We found 
1/4 LRGs physically associated with the Mgll absorber 
and all four pairs are within p = 100 kpc . We are 
currently increasing the number of spectroscopic pairs 
from Sample B to obtain a more accurate measurement 
of K. 

5. DISCUSSION 

Our follow-up spectroscopic study of close LRG-Mg II 
pairs has confirmed a physical association between the 
Mgll absorber and the LRG in five of the 15 cases 
studied. The small velocity separations, \Av\ < 350 
kms~^ and projected distances p < 160 kpc demon- 
strate that at least some of the strong Mgll absorbers 
[Wr(2796) > 1 A] that contributed to the small-scale 
[rp = {1 -\- z) X p < 600 kpc] clustering signal in 
Gauthier et al. (2009) do indeed originate in the hosting 
dark matter halos of LRGs. According to a stell ar popu- 
lation analysis ( using the stellar templates from Bruzual 
fc Charlot|2003 ) done on the LRG spectra, we found that 
these galaxies are characterized by a stellar population 
at least IGyr old with a most likely age between 1-11.75 
Gyr. The optical spectra are characteristic of an evolved 
stellar population with little recent star formation. 

In addition, our search of Mg II absorption features in 
the vicinity of LRGs has yielded a conservative estimate 
of the maximum covering fraction of cool gas in LRG ha- 
los. Because of the redshift uncertainties of the LRGs, we 
cannot distinguish between whether the observed Mgll 
absorbers are physically associated with the LRG ha- 
los or these absorbers occur due to chance coincidence, 
correlated galaxies, or surrounding satellite galaxies. De- 
spite these various uncertainties, the observed incidence 
of Mg II absorbers places a strong constraint on the max- 
imum covering fraction of cool gas in massive LRG halos 
as probed by the Mgll absorbers. We find that the cool 
gas covering fraction of absorbers with VF^(2796) > 0.5 
A in massive LRG halos is no more than 20% and no 
more than 7% for stronger absorbers VF^(2796) > 1 A. 

Although to establish a physical association between 
these close LRG-Mg II absorber pairs requires spectro- 
scopic data of the LRGs, we assess possible contami- 
nations due to surrounding satellite galaxies and back- 
ground random/correlated absorbers based on theoreti- 
cal expectations of the satellite population and the mean 
number density of Mgll absorbers per line of sight ob- 
served along random QSO sight lines. 

To estimate the incidence of Mgll absorbers due to 
satellite galaxies in the LRG halos, we first adopt the 
mass function of subhalos (Equation 10 in Tinker & Wet-] 
zel 2009). The mass of subhalos is defined as the mass 
at the time of accretion — this mass will correlate with 
the gaseous halo that the subhalo initially had. Next, we 
adopt the gas radius at Wr{2796) > 1 A and the covering 
fract ion within the gas radiu s as a function of halo mass 
from [Tinker & Chen ( |2010 |). The total cross section of 
subhalos within host halos of mass 10^^ h-^ Mq is then 
estimated according to 



^sub 



1 



O-(Mhost) 

XTrR\M,^k,Wr 
where Mhost is the host halo mass (in this case, 10^^ 



cZM,ubn(M,ub|Mhost) (1) 

l)Kg(Msub) 



Incidence of cool gas in 10^^ halos 



7 



h ^ Mq ) , Msub is subhalo mass at time of ac cretion, 
n(Msubl-^host) is the subhalo mass function from [Tinker 



fc Wetzel (2009), R is the gas radius for Wr {27m) 1 



A and K.g is the covering fractio n within the gas radius . 
Both R and K/g are adopted from [Tinker fc Chen| ( [2QTq| , 
which are calibrated on data to match the frequency 
function of absorbers and large scale bias as a function 
of Vl^^ (2796). 

Under the assumption that the gaseous halos of satel- 
lite halos are the same as field halos of similar mass, 
Equation 1 yields a maximal covering fraction of 1.8% 
at p < 350 kpc . This is the maximum value be- 
cause tidal str ipping and ram pressure are n ot taken into 
account (e.g., [d unn fc Got t[[1972[ [Balogh et al. 2000). 
Equation 1 is also maximal because it assumes that the 
subhalos have a random distribution within the parent 
halo. If we further assume that the subhalos follow an 
NFW profile (Wetzel & White 2009), ^sub depends on 
impact parameter as shown in Figure 4. The satellite 
contribution is at most comparable to the A^max mea- 
surements at p < 100 kpc , but the satellite cover- 
ing factor of cool gas is > 1 dex below /^max at p > 200 
kpc . 

To estimate the incidence of Mgll absorbers due to 
structure along a line of sight that intersects a large 
halo, we adopt the observed mean number density of 
W^(2796) > 1 A absorbers from Nestor et al. (2005) and 
Prochter et al. (2006). These authors report a number 
density of n ^ 0.1 per unit redshift interval per line of 
sight at z = 0.5 for Mgll absorbers of W^(2796) > 1 A. 
Over a redshift interval of Az = 0.27, which corresponds 
to the redshift uncertainty of the LRGs, we therefore 
expect a background contamination of 0.027 absorbers 
per line of sight. To estimate the incidence of Mgll ab- 
sorbers due to correlated large-scale galaxy distribution, 
we calculate the predicted incidence of correlated ha- 
los around the LRGs. We select halos becau se they 



are e xpected to be the dominant contaminant (Tinker 
fc Ch en 20081 [2 010[). First, we adopt Equation (11) in 



TmEer fc Che^J(^ 



0) and modify this equation to take 

into account the correlated structures along the line of 
sight. We integrate along the line-of-sight distance from 
—300 to +300 Mpc . We repeat the calculation for 
each projected distance to the LRG and determine the 
number of halos that could contribute to the signal based 
on their mass and size. Our calculation indicates roughly 
10% enhancement in the background number density of 
Mgll absorbers per line of sight. We therefore conclude 
that background random and correlated galaxies would 
have contributed no more than 3% of the observed Mg II 
absorbers around the massive LRGs. 

Our exercise shows that contaminations due to 



surrounding satellite galaxies and background ran- 
dom/correlated galaxies are insufficient to account for 
the observed incidence of Mgll absorbers in the vicin- 
ity of LRGs. The maximum covering fraction of strong 
absorbers is ~ 7% averaged over the entire halo. 

fBowen et al. (2006) reported 100% coverine; fraction of 
Mgll absorbers within 100 kpc projected distance from 
four foreground QSO at z = 0.65 — 1.55. The observed 
large incidence of Mgll absorbers in QSO halos appears 
to be inconsistent with the low gas covering fraction 
found in the LRG halos. This discrepancy can be un- 
derstood b y considering the mean mass scales of QSOs 
and LRGs. Ross et al. (2009) calculated the clustering of 
SDSS spectroscopic QSOs at z < 2.2. They found that, 
at z ~ 0.5, the bias of the dark matter halos traced by 
QSOs is 6 ~ 1.4 which is smaller than th e bias estimate 
for th e LRGs at the same redshift {b ~ 2; [Gauthier et al 
20091). This implies that the typical mass of the dark 



matter halos hos ting QSOs is aro und log M/h ~ 12.3 
(see figure 12 in Ross et al. 2009), which is about one 
dex smaller than what is founS for LRGs in IGauthierl 



et al.j ( 2009 ) . The difference in the observed gas covering 
fraction therefore reflects the differences in the c old gas 



content of halos on different mass scales (see e.g., [Tinker 
'fc Chen'2008, Chen et al. 2010 - ApJ submitted). 

Recall that Mgll absorbers trace high column density 
gas of neutral hydroge n column density A/' (HI) > 10^^ 
cm~^ (Bergeron < 



Stasihska|l986|[Rao et al.|2006[ ). Fur 
l1. (2006) nave shown that on average 



thermore, Rao et al. (2006) nave snown tnat on average 
42% of iy^(2796) > 0.6 A Mgll absorbers with asso- 
ciated strong Fell and Mgl absorption transitions are 
likely to contain neutral gas of A^(HI) > 2 x 10^^ cm~^. 
Two of the Mgll absorbers presented in Figure 2 sat- 
isfy these criteria and are therefore likely damped Lya 
absorbers. Our study therefore provide direct empirical 
evidence for the presence of high column density cool gas 
in 10^^ Mq halos. A large spectroscopic sample of 
LRGs will not only improve the contraints for Kg in mas- 
sive halos but also provide additional insights into the 
origin of the cool gas. 
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